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We explore three-nucleon forces (3NF) from lattice QCD simulations. Utilizing the Nambu- 
Bethe-Salpeter (NBS) wave function, two-nucleon forces (2NF) and 3NF are determined on the 
same footing. Quantum numbers of the three-nucleon (3N) system are chosen to be (/,/^) = 
(1/2,1/2+) (the triton channel). The enormous computational cost is reduced by employing 
the simplest geometrical configuration, where 3N are aligned linearly with an equal spacing. We 
perform lattice QCD simulations using Nf = 2 dynamical clover fermion configurations generated 
by CP-PACS Collaboration, at the lattice spacing of a = 0. 156 fm on a 16^ x 32 lattice with a large 
quark mass corresponding to nij^ = 1.13 GeV. Repulsive 3NF is found at short distance. 
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1. Introduction 



Determination of the properties of three-nucleon forces (3NF) is one of the most important 
issue in nuclear physics and astrophysics these days. There are various phenomena where 3NF are 
considered to play an important role, e.g., the binding energies of light nuclei [jl|], the properties of 
neutron-rich nuclei and the supernova nucleosynthesis [^, the nuclear equation of state (EoS) at 
high density relevant to the physics of neutron stars [Q, ^ and the cross sections in nucleus-nucleus 
elastic scatterings Deuteron-proton elastic scattering experiments at intermediate energies have 
also shown a clear indication of 3NF 

Despite of its phenomenological importance, microscopic understanding of 3NF is still lim- 
ited. Pioneered by Fujita and Miyazawa the long range part of 3NF has been commonly mod- 
eled by the two-pion exchange (27rE), particularly with the A-resonance excitation. An additional 
repulsive component of 3NF at short distance is often introduced in a purely phenomenological 
way [^. An approach based on the chiral effective field theory is also developing [^. 

However, since 3NF is originated by the fact that a nucleon is not a fundamental particle, it 
is most desirable to determine 3NF directly from the fundamental degrees of freedom (DoF), i.e., 
quarks and gluons, on the basis of quantum chromodynamics (QCD). In this proceeding, we carry 
out first-principle calculations of 3NF using lattice QCD simulations. Note that while there are 
lattice QCD works for three- and four- baryon systems [ p^ , 11], they focus on the energies of the 
multi-baryon systems, and extracting 3NF is currently beyond their scope. 

As for the calculation of two-nucleon forces (2NF) from lattice QCD, an approach based on 

the NBS wave function has been proposed [12, 13], so that the potential is faithful to the phase 

shift by construction. Resultant (parity-even) 2NF are found to have desirable features such as 

attractive wells at long and medium distances and central repulsive cores at short distance. The 

method has been successfully extended to the hyperon-nucleon (YN) and hyperon-hyperon (YY) 

interactions [|I4KI|, |T6|, 1^, [iSl] and meson-baryon systems []19|]. In this report, we extend the 

method to the three-nucleon (3N) system, and perform the lattice QCD simulations for 3NF in the 
jP 



triton channel, (7,7^) = (1/2, 1 /2+) [Eg, 11, 22, M]. For details of this study, refer to Ref. [22]. 



2. Formalism 

The detailed formulation for the determination of 2NF is given in Ref. [[l3|], and we here focus 
on the extension to the 3N system. We consider the NBS wave function of the 3N, Y3N{r,p), 
extracted from the six-point correlator as 

G3Nir,p,t-to) = ^ X^(0| (A^(xi )N{x2)Nix,)){t) {N'N'N') {to)\0) (2.1) 

R 

> A^N¥3Nir,p)e-^"'^'-'°\ (2.2) 

f>;o 

\lf3Nir,p) = {0\Nixy)N{x2)N{x3)\E3N), A3n^{E3n\{N'N'N')\0), (2.3) 

where Et,n and l^sAr) denote the energy and the state vector of the 3N ground state, respectively, N 
(N') the nucleon operator in the sink (source), and R= {x\ +X2+XT,)/3,r = x\ —X2, p =X3 — {xi + 
X2)/2 the Jacobi coordinates. 
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With the derivative expansion of the potentials |24], the NBS wave function can be converted 
to the potentials through the following Schrodinger equation, 

1 1 



2Hp 



yl+Y.^2Ninj) + V3NF{7,p) 



(2.4) 



where VjNirij) with r,y = Xj — xj denotes the 2NF between (i,7)-pair, VsA^f (r,p) the 3NF, /i^ = 
niN/l, Hp = 2mj^/2> the reduced masses. If we calculate ^f-ii^{f,p) for all rand p, and if all V2w(r,y) 
are obtained by (separate) lattice calculations for genuine 2N systems, we can extract Vj,i^f{f,p) 
through Eq. ( ^ . 

In practice, however, the computational cost is enormous, because of enlarged color/spinor 
DoF by the 3N (i.e., 9 valence quarks) and factorial number of the Wick contractions. In order 
to reduce the cost, we first take advantage of symmetries (such as isospin symmetry) to reduce 
the number of the Wick contractions. Second, we utilize that there is a freedom for the choice 
of a nucleon interpolating operator. In particular, a potential is independent of the choice of a 
nucleon operator at the source, N' , and we employ the non-relativistic limit operator as A^' = Nnr = 
£ahc{(laCY5P„r(]h)Pnr(lc with = (1 + Ya) /2. Compared to the standard nucleon operator, Ngtd = 
£abc{4^Cysqb)qc, the spinor DoF of each (source) nucleon is reduced by half, and the computational 
cost of the 3N system is reduced by a factor of 2^ = 8. On the other hand, a potential is dependent 
on the choice of a nucleon operator at the sink, N, since a NBS wave function is defined through a 
sink operator. Note, however, that physical observables calculated from these different potentials, 
such as phase shifts and binding energies, are unique by construction [13]. In this sense, choosing 
corresponds to choosing the "scheme" to define the potential [13, 25]. In Ref. [^], it is found that 
the non-locality of 2NF from the choice ofN= N^td is small, and N = Nstd can be considered to be 
a "good scheme". We therefore employ the same sink operator A'^ = Nstd in the 3N study as well, so 
that 2NF and 3NF are determined on the same footing. Finally, we restrict the geometry of the 3N. 
More specifically, we consider the "linear setup"with p = 0, with which 3N are aligned linearly 
with equal spacings of r2 = \r\/2. In this setup, the third nucleon is attached to (l,2)-nucleon 
pair with only S-wave. Considering the total 3N quantum numbers of (/,/^) = (1/2,1/2+), the 
triton channel, the wave function can be completely spanned by only three bases, which can be 
labeled by the quantum numbers of (l,2)-pair as '^o, ^^i, ^D\. Therefore, the Schrodinger equation 
leads to the 3x3 coupled channel equations with the bases of ^\s^, Y^sr V^'Dr The reduction 
of the dimension of bases is expected to improve the S/N as well. It is worth mentioning that 
considering the linear setup is not an approximation: Among various geometric components of 
the wave function of the ground state in the triton channel, we calculate the (exact) linear setup 
component as a convenient choice to study 3NF. While we can access only a part of 3NF from it, 
we plan to extend the calculation to more general geometries step by step, toward the complete 
determination of the full 3NF. 

We consider the identification of genuine 3NF. It is a nontrivial work: Although both of parity- 
even and parity-odd 2NF are required to subtract 2NF part in Eq. (2.4), parity-odd 2NF have not 
been obtained yet in lattice QCD. (See, however, our recent progress on this issue [26].) In order 
to resolve this issue, we consider the following channel, 

1 



^^5 = 



(2.5) 



3 



Three-Nucleon Forces explored by Lattice QCD Simulations 



Takumi Doi 



which is anti-symmetric in spin/isospin spaces for any 2N-pair. Combined with the Pauli-principle, 
it is automatically guaranteed that any 2N-pair couples with even parity only. Therefore, we can 
extract 3NF unambiguously using only parity-even 2NF. Note that no assumption on the choice 
of SD-configuration of r, p is imposed in this argument, and we thus can take advantage of this 
feature for future 3NF calculations with various setup of 3D-geometries. 



3. Lattice QCD setup and Numerical results 

We employ Nf = 2 dynamical configurations with mean field improved clover fermion and 
RG-improved gauge action generated by CP-PACS Collaboration [p7|]. We use 598 configurations 
at j3 = 1.95 and the lattice spacing of a^^ = 1.269(14) GeV, and the lattice size ofV =L^ xT = 
16^ X 32 corresponds to (2.5 fm)^ box in physical spacial size. For u, d quark masses, we take the 
hopping parameter at the unitary point as Kud = 0.13750, which corresponds to mj^ = 1.13 GeV, 
m^? = 2.15 GeV and = 2.31 GeV. We use the wall quark source with Coulomb gauge fixing. 
In order to enhance the statistics, we perform the measurement on 32 wall sources using different 
time slices, and the forward and backward propagations are averaged. The results from both of 
total angular momentum 7, = ±1/2 are averaged as well. Due to the enormous computational cost, 
we can perform the simulations only at a few sink time slices. Looking for the range of sink time 
where the ground state saturation is achieved, we carry out preparatory simulations for effective 
2NF in the 3N system in the triton channel at 2 < (? — ?o)/a < 11, and find that the results are 
consistent with each other as long as {t — to) /a > 1 [^]. Being on the safer side, we perform linear 
setup calculations of 3NF at {t — to)/a = 8 and 9. We perform the simulation at eleven physical 
points of the distance r2. 

In Fig. [l|, we plot the radial part of each wave function of Ys = {—V'So ~^ V'^sJ/V^, 1/% = 
iV'So~^ i//35j)/\/2 and y^Di obtained at {t — Iq) / a = 8. Here, we normalize the wave functions by 
the central value of Ysi^i = 0)- What is noteworthy is that the wave functions are obtained with 
good precision, which is quite nontrivial for the 3N system. We observe that y/^ overwhelms other 
wave functions. This indicates that higher partial wave components are strongly suppressed, and 
thus the effect of the next leading order in the derivative expansion, spin-orbit forces, is suppressed 
in this lattice setup. 

We determine 3NF by subtracting 2NF from total potentials in the 3N system. Since we have 



only one channel (Eq. (|2.5[)) which is free from parity-odd 2NF, we can determine one type of 3NF. 
In this report, 3NF are effectively represented in a scalar-isoscalar functional form, which is often 
employed for the short-range 3NF in phenomenology [^]. 

In Fig. ^, we plot the results for the effective scalar-isoscalar 3NF at {t — to)/a = 8. Here, we 
include r2 -independent shift by energies, 5e — 5 MeV, which is determined by long-range behavior 



of potentials (2NF and effective 2NF in the 3N system) While Se suffers from ^10 MeV 
systematic error, it does not affect the following discussions much, since 5e merely serves as an 
overall offset. In order to check the dependence on the sink time slice, we compare 3NF from 
{t — to)/a = ^ and 9 in Fig. |3[ While the results with {t — to)/a = 9 suffer from quite large errors, 
they are consistent with each other within statistical fluctuations. 

Fig. ^shows that 3NF are small at the long distance region of ri- This is in accordance with the 
suppression of 27rE-3NF by the heavy pion. At the short distance region, however, an indication of 
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Figure 1: 3N wave functions at (f — to) /a — 8. 
Circle (red), triangle (blue), square (green) points 
denote y/s, V^Di ' respectively. 




Figure 2: The effective scalar-isoscalar 3NF in 
the triton channel with the linear setup, obtained at 
(f-fo)/fl = 8. 



repulsive 3NF is observed. Note that a repulsive short-range 3NF is phenomenologically required to 
explain the properties of high density matter. Since multi-meson exchanges are strongly suppressed 
by the large quark mass, the origin of this short-range 3NF may be attributed to the quark and 
gluon dynamics directly. In fact, we recall that the short-range repulsive (or attractive) cores in the 
generalized two-baryon potentials are systematically calculated in lattice QCD in the flavor SU(3) 
limit, and the results are found to be well explained from the viewpoint of the Pauli exclusion 
principle in the quark level [|l^]. In this context, it is intuitive to expect that the 3N system is subject 
to extra Pauli repulsion effect, which could be an origin of the observed short-range repulsive 3NF. 
Further investigation along this line is certainly an interesting subject in future. 

It is in order to discuss the systematic errors in lattice simulations. First, one may worry about 
the discretization error, since the nontrivial results are obtained at short distance. In particular, the 
kinetic term (Laplacian part in the Schrodinger equation) could suffer from a substantial effect, 
since they are calculated by the finite difference as V^jj/(x) = ^ [f{x + a,) + f{x — Ui)— 2/(x)] . 
In order to estimate this artifact, we also analyze using the improved Laplacian operator for both of 
2N and 3N, Wl^^f{x) = I,- [- (/(x + la^) + /(x - 2a,-) ) + 16(/(x + a,-) + f{x - a,-) ) - 30/(x) )] . 
In Fig. ^ we plot the comparison between the results of 3NF from V^^^ and V?^p at {t — to)/a = 8. 
It is found that they are consistent with each other, and we conclude that the discretization artifact 
of 3NF in Laplacian operator is small. We, however, remark that this study probes only a part of 
discretization errors. Actually, the analysis with operator product expansion ^ shows that 
2NF of V2N{r) tend to diverge as r — 0, so significant discretization artifact is expected around 
r = 0. Full account of the discretization artifact can be examined by an explicit lattice simulation 
with a finer lattice, which is currently underway. 

Second, the finite volume artifact is discussed. In this simulation, three nucleons are accom- 
modated in (2.5 fm)-' spacial lattice box. We note that this is quite a large box for the heavy pion, 
namely, nijcL is as large as 14. Furthermore, the finite volume artifact is expected to appear mainly 
in large r2 region, and we avoid it as much as possible by focusing on the short-range part of 3NF. 
For the relatively large r2 (r2 ^ 0.5 fm) region, points are carefully chosen so that they are located 
in off-axis directions. Therefore, we expect that the finite volume artifact is not substantial in our 
study. Of course, a quantitative examination requires calculations with larger volumes, which we 
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Figure 4: 3NF obtained at (f - fo)/a = 8. Cir- 
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defer to future studies. 

Third, we consider the contamination from excited states. As has been discussed, we do 
not observe sink time dependence for 3NF in hnear setup between {t — toj/a = ^ and 9, nor for 
effective 2NF in the 3N system as long as {t — to)/a > 7. It is, however, desirable to investigate 



3NF in linear setup with more sink time slices. In particular, it is recently proposed [ [30| ] to utilize 
the time-dependent Schrodinger equation to further eliminate the excited state contamination. In 
order to apply this method to 3NF, linear setup calculations at additional sink time slices are in 
progress. 

Finally, quark mass dependence of 3NF is certainly an important issue, since the lattice simula- 
tions are carried out only at single large quark mass. In the case of 2NF, short-range cores have the 



enhanced strength and broaden range by decreasing the quark mass []130. We, therefore, would ex- 
pect a significant quark mass dependence exist in short-range 3NF as well. In addition, long-range 
2;rE-3NF will emerge at lighter quark masses, in particular, at the physical point. Quantitative 
investigation through lattice simulations with lighter quark masses are currently underway. 
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